Investigation of vibration’s effect on driver in optimal motion cueing algorithm

The increased sensation error between the surroundings and the driver is a major problem in driving simulators, resulting in unrealistic motion cues. Intelligent control schemes have to be developed to provide realistic motion cues to the driver. The driver’s body model incorporates the effects of vibrations on the driver’s health, comfort, perception, and motion sickness, and most of the current research on motion cueing has not considered these factors. This article proposes a novel optimal motion cueing algorithm that utilizes the driver’s body model in conjunction with the driver’s perception model to minimize the sensation error. Moreover, this article employs H∞ control in place of the linear quadratic regulator to optimize the quadratic cost function of sensation error. As compared to state of the art, we achieve decreased sensation error in terms of small root-mean-square difference (70%, 61%, and 84% decrease in case of longitudinal acceleration, lateral acceleration, and yaw velocity, respectively) and improved coefficient of cross-correlation (3% and 1% increase in case of longitudinal and lateral acceleration, respectively).


Introduction
Evaluation of Virtual Reality in driving and on-road field studies have shown encouraging resemblances [1,2].However, when compared to natural world driving, questions are raised on the validity of driving simulators.It will not be easy to perceive a driver's behavior in such an environment.Driving simulators are considered very important to improve road safety through proper indoor training and improved vehicle design [3].Their use is crucial when the number of vehicles is considerably large, and road structure is also not very appropriate.Research has shown that a virtually trained driver saves money, time, life, and material.Nowadays, driving simulators are used not only as training simulators but also as entertainment simulators.Researchers have evaluated drivers' behavior, performance, and attention.In the automobile industry, they are used for optimal and risk-free vehicle design, resulting in reduced vehicle cost and increased safety.We can have efficient hardware implementations of driving simulators using energy efficiency techniques, such as SuperSlash [4], which is used for energy efficiency of artificial intelligence (AI) based models.
The motion cueing algorithm (MCA) converts linear accelerations and angular velocities of a virtual vehicle into movements of a driving simulator.One of the first methods to transform these accelerations and velocities was the classical MCA proposed by Conrad et al. [5].Classical MCA is a linear cueing algorithm that is simple and easily tunable and has been widely used in different types of simulators [1,6].Classical motion cueing has a short processing time and stable performance [7].However, classical MCA has many drawbacks.For example, it needs a regulation procedure that emphasizes the worst-case situation resulting in conventional motion with deprived workspace usage, hence not being flexible in all situations.Moreover, it completely ignores the human vestibular system [8] and may generate false cues, distort the signal, and lead to wrong motion signals due to filter characteristics [9].Also, a lowpass filter leads to phase delay [10].The coordinated adaptive washout algorithm has been applied to flight simulators in [11].The authors replaced the tilt coordination block of the classical washout filter with an adaptive filter to avoid residual motion base displacements.In [12], adaptive motion cueing using fuzzy-based tilt coordination has been proposed.Due to many other possibilities of tilt coordination [13], the study of this algorithm is just at its initial stages.
An optimal MCA was proposed ignoring the vestibular system [14].Optimal MCA has some drawbacks, but it generally reduces the sensation error between the virtual surrounding and real drivers.However, the results are not satisfactory to a real human perception due to physical limitations and an optimal washout filter needs to be applied to a larger workspace [15].Getting realistic motion cues on a six-degree-of-freedom Stewart platform compared to the actual vehicle is still a problem.In some cases, the motion system has to be effectively turned off to avoid improper motion cues.A robust optimal MCA based on the linear quadratic regulator (LQR) method and the genetic algorithm (GA) has been studied in [16].The driver feedback has been used in [17] to characterize desirable features of race cars; then the numerical optimal control (NOC) was used to decide platform motions.NOC is a useful open-loop method for analyzing performance constraints of the driver in loop strategies [17].Other methods, including LQR-based optimal control and GA-based optimal control [16], exhibit better performance in improving the sensation output, producing realistic motion cues but cannot enhance workspace usage due to several constraints.The LQR-based method improves the workspace, but it has more sensation error and the GA-based method decreases error but with poor workspace usage [16].Recently, many methods have introduced neural networks and computational intelligence-based motion cueing [18,19]; however, these techniques are computationally complex.
This research utilizes H1-based optimal MCA, which shows superior performance compared to the state-of-the-art methods.This article has the following main contributions: 1.A novel optimal MCA is proposed that utilizes the driver's body model together with the driver's perception model.
2. The LQR is replaced with an H1-based controller to optimize the quadratic cost function of sensation error.
A significant advantage of using H1-based MCA is that it can be used for any driving simulator without parameter tuning.We need to adjust the physical constraints of the new system and do not need to care about "trial and error" for parameter tuning.The article is organized as follows: the Background section comprises the math behind the Stewart platform and its kinematics; the Proposed Methodology section includes the application of proposed optimal MCA on the Stewart platform; the Results and Discussion section compares different motion cueing techniques.Finally, the Conclusion section concludes the article by discussing the improvements made by the proposed methodology.

Background
In this section, we will discuss the mathematics behind the Stewart platform along with its kinematics.Fig 1 refers to the schematics of a generic Stewart platform.There are two parts to every limb length (l i ) connected by a prismatic joint (P i ).This limb's upper and lower sides are connected to the upper and lower platforms with spherical (or universal) joints.In this robot, six identical hydraulic actuators can be actuated to obtain six-degree-of-freedom motions of the moving platform.Six universal joints at the base platform of the Stewart platform are referred to as A i , where i = 1, 2, . . .6. Similarly, six spherical joints at the moving platform of the Stewart platform are referred to as B i , where i = 1, 2, . . .6.The prismatic joints between fixed and moving platforms are referred to as P i , where i = 1, 2, . . .6.The limb length of every kinematic chain of Stewart platform is referred to as li , where i = 1, 2, . . .6. Two frames O A and O B are attached to the fixed and moving platforms, respectively.

Inverse kinematics
Inverse kinematics deals with the determination of limb lengths when the position and orientation of the moving platform are known [20][21][22].In the case of a parallel robot, like a Stewart platform, inverse kinematics is easy [23].The position AP of the moving platform and orientation ŝ are given, and limb lengths l i s are to be determined.The position is given as follows: In the case of screw axis representation, the orientation of the platform is given as follows: In the case of rotation matrix representation, the orientation of the platform is given as follows: Using the loop closure equation for each limb on the geometry of Fig 1, The procedure in [17] will be followed to find actuator velocity.Taking the derivative of Eq (5) and simplifying for angular velocity, we get

Forward kinematics
Forward kinematics is the method of determining the position and orientation of moving platforms when limb lengths of the Stewart platform are known.As shown in Fig 1, if vector L of limb lengths l i , where i = 1,2,� � �,6, is known, then the problem is to determine position A P and orientation A B R. The optimal solution method having many iterations needs high-speed computation due to the heavy computational burden [24].A good initial guess has to be made to establish a rapid and converging solution.Forward kinematics can be obtained using nonlinear least-squares algorithms based on Levenberg-Marquardt or trust-region-reflective methods [23].The trustregion-reflective method is based on the interior-reflective Newton method proposed in [25].The Levenberg-Marquardt method has been proposed in [26,27].To apply one of these algorithms, the first step will be to merge Eqs (1) and (2): or, equivalently, Position A P and orientation A B R will be determined as follows: where vt 7 = (1-cos (t 7 )), st 7 = sin (t 7 ) and ct 7 = cos (t7).
From Fig 1, the moving platform has points B i 's and the base platform has points A i 's and their vectors with respect to the base are b i 's and a i 's, respectively.The following set of nonlinear equations can be solved for the forward kinematics: and Optimization methods can be employed here to solve the above seven nonlinear equations (Eqs ( 12) and ( 13)).In case of least-squares optimization, minimize f ðtÞ ¼ The solution of forward kinematics is demonstrated as a flowchart in Fig 2 .In the first flowchart process, Eqs 9 and 10 are used for the calculation of p and R, respectively.Then, block F i is evaluated using Eq (12).Now, f(t) is determined as demonstrated above.The conditional block checks if this value is less than required accuracy ε or not.In case f(t) is less than ε, the function iteratively calculates the new value of t and sends it back to calculate p and R, to recalculate the new values of p and R corresponding to the new value of t.
The input trajectories of position and orientation are shown in

Motion cueing algorithm
Motion cueing refers to the approaches through which the environment of a driving simulator is made to feel like a real environment.These approaches are known as cues and are classified as acoustic, visual, and motion cues.A driving simulator of any class must contain an efficient visual system.Moreover, most visual systems are accompanied by acoustic systems.Motion simulators are required to provide motion cues sensed by vestibular, body, and perception system models.In the real environment, motion system parameter values exceed the limits of the driving simulator motion system.Therefore, it is essential to convert those parameters.The methods applied to perform this conversion are known as MCAs, motion drive algorithms, or washout filters.These algorithms aim to reproduce the motion cues in a realistic way, ensuring the usage of maximum workspace.The first was classical motion cueing built in [5] for flight simulators.Naturally, the requirements for a driving simulator are somewhat different; nonetheless, the classical algorithm can be adapted with some modifications.

Proposed methodology
Optimal MCA for vehicle driving simulators that incorporate driver's body and perception models is proposed so as to include human body vibrations [15,[28][29][30].Human vibrations are of three types: the mass-spring-damper model [31], the model from the finite element method [32], and the transfer function model [14].Methods for measuring periodic, random, and transient whole-body vibrations have been discussed in [28].
Fig 5 shows the basicentric coordinate system of a human body commonly used to illustrate human body vibrations.The driver's body model is given by Eqs (14-19) as a transfer function that contains a band-limiting filter GBL in Eq (18).This filter contains high-pass filter G HP where a i , i = 1, 2, 3, . . .,6 are the frequency weighted acceleration values and H i , i = 1, 2, 3, . . .,6 are the resonant quality factors.These factors evaluate the overall frequency weighting of the transfer function.The values of these variables are obtained from [28].The driver's perception model comprises vestibular and otolith organs.These organs are responsible for the perception of rotational velocity and specific force.Specific force is the difference between linear acceleration and gravity.Two transfer functions are shown in Eqs (20) and (21).The first transfer function is the ratio of the vestibular organ to that of rotational velocity: The second transfer function is the ratio of the otolith organ to the specific force given as follows: where L L , L A , and L S are the weighting factors from vestibular organ and τ L , τ A , and τ S are weighting factors from the otolith organ.

Problem statement
The block diagram of a driving simulator, including an optimal motion cueing block W(s), is shown in Fig 6.
The input to this block is u v (s) and its output is u s (s), which serves as input to the driver's body model.The six motions of the moving platform are grouped into four modes as follows: 1. Longitudinal mode 2. Lateral mode 3. Heave mode

Yaw mode
The longitudinal mode contains surge and pitch as linear and angular motion components, respectively.Sway has linear motion and corresponding roll as an angular motion component in the case of the lateral mode.Heave mode only contains heave as a linear motion component and yaw is the only component in yaw mode.In each mode, the state-space equation can be obtained from the transfer function of the driver's body models.The following text includes the description of the algorithm of longitudinal mode, where Ω p is the angular velocity or derivative of pitch and a x is linear acceleration along with surge.
The structure of the proposed MCA is given in Fig 6 .Eq (22) shows longitudinal components of input u v (s) to simulator block.
The specific force f x (s) is defined as force per unit mass along surge and is given as follows: The output of this cueing model is u s (s): The problem is to determine a transfer function W s (s) that will be used to compute the output sensational signals for simulator moving platform position and orientation.
Determination of sensation signal.In optimal control theory, a cost function is to be optimized considering the motion system constraints [32][33][34].In this case, the sensation error between virtual surroundings and real drivers is minimized.Sensation error can be determined by the observable canonical form of state space.In the case of angular motion channel, Eqs (25) and ( 26) are used, and in the case of translational motion channel, Eqs ( 27) and ( 28) are used.
Similarly, using Eq (23), Combining Eqs ( 25)- (28) with Eqs ( 29) and ( 30), If the input to driver's body model is u s , then similar equations are obtained: Error dynamic equation can be obtained from Eqs ( 29)- (32) and is given as follows: It is seen in practical situations that the vehicle signal is accompanied by a filtered white noise signal, so it must be considered in the error signal.Eqs (35) and (36) show the state space of vehicle signal for noise considering w as filtered white noise input: Adding the effect of white noise to the error dynamic equation for the determination of cumulative error e, where The main objective of this article is to produce realistic motion cues, i.e., motion cues similar to the real vehicle.To achieve this objective, the error between the vehicle signal and simulator motion signal, given by Eq (38), must be minimized.In the last section error, the dynamic equation is modeled; in this section, this error is minimized using the appropriate choice of the optimization method.As explained in the motion cueing section, there are workspace limitations of the motion simulator.These limitations must be included in the optimal control problem to generate realistic motion cues.Objective and constraints of motion system can be combined into a quadratic cost function and the problem is formulated as follows: Find value of u s that minimizes J: Subject to the following constraints, where J is the quadratic cost function to be minimized, Q and R are weighting matrices, and l i ,v i ,a i ,i = 1,2,3,� � �6 are limb length, velocity, and acceleration of every leg of the motion system, respectively.
The enhanced sensitivity to the sensation error.The problem is to decide the inputs u s such that our quadratic cost function is minimized within the limits of constraints.Fig 6 shows the input vehicle dynamics u v as a vector of longitudinal acceleration, lateral acceleration, and yaw velocity.A similar vector u s is generated as simulator input.Both of these inputs are used for the calculation of states of the system and outputs as discussed in Eq ( 30 In short, there is a need of selecting such parameters which minimize the sensed longitudinal acceleration, lateral acceleration, and yaw velocity.The accelerations are translational whereas the velocity is a rotational parameter.However, this is a lengthy iterative process and requires repeated driver feedback.This work includes driver's body model, to minimize the workload of drivers as shown in Eq (46) below: Analysis of cost function.We require to cue the brake system and the accelerate system while in longitudinal mode.For both situations we use numerical control techniques for error free acceleration profiling.The basicentric coordinate system of a human body in the maneuver at the front of the workspace and must conclude it there for the optimum reaction to operate.As a result, the platform was set up for braking, and the maneuver concluded in an acceleration-ready posture.This cueing signal has two characteristics that need discussion.The first is a delayed reaction that causes the platform's peak deceleration to coincide with the cars.The platform comes to a stop near the end of the braking action due to a mistake made during workspace management, which is the motion's second distinguishing feature.This answer falls short of what the racing car drivers we contacted demand in terms of a strong early onset cue, hence it is not acceptable.Discussions led to the conclusion that drivers need a strong first onset cue followed by a mistake of the smallest feasible size.Since the cueing acceleration produced lacks these traits, it needed to be modified.
In a drawn-out high-speed turn, lateral acceleration signals can be of great amplitude and extended duration and have a substantially different form from the longitudinal case.Simulators with a constrained workspace cannot simulate this kind of acceleration.Again, the acceleration at the beginning is cued but not maintained.Since the yaw acceleration and the lateral cue have a similar shape, an objective function of the type Eq (44) is also utilized in this situation.

Results and discussion
A good motion cueing system should possess the following properties: 1. Realistic motion cues, i.e., the wave shape of output acceleration cues, should be more similar to that of input linear acceleration.Similarly, the wave shape of output angular velocity cues should be more similar to that of input angular velocity.
2. Usage of more workspace, i.e., the limits of theoretical workspace, should be achieved by optimum position and orientation of motion base.In this article, all the above three qualities have been practically achieved.Results are compared based on the root-mean-square difference (RMSD) and coefficient of cross-correlation (COCC) that show similarity between two waveforms to produce realistic motion cues.The following section includes the application of control strategies step by step for motion cueing and compares the results.

State-of-the-art motion cueing algorithms
Optimal MCA based on the linear quadratic regulator (LQR) is modified by changing the axis of rotation to the driver's body axis instead of the motion system's center of rotation [16,22].The main aim behind this change is to achieve realistic motion cues and remove false cues.Fig 5 shows the basicentric coordinate system of the driver body (seated position).The cross-coupling between the simulator and driver depends on the distance between their respective coordinate systems.As discussed in the motion cueing section, MCA receives linear acceleration and angular velocity as input to produce a proper cueing signal for the simulator motion base.Furthermore, a tilt rate limiter is introduced after cueing has been achieved for better tilt coordination.This addition investigates the effects on both motion cue fidelity and displacement.
The translational input has time values of longitudinal, lateral, and heave-specific forces.These inputs vary depending on acceleration, deceleration, and retardation, as shown in Figs 7  and 8.In these Figs., the blue plot shows values of input linear acceleration and illustrates that these values are never uniform; i.e., they keep on changing with time.According to Fig 7, longitudinal specific force values are applied to linear quadratic regulator-based optimal MCA, and resultant cues are plotted in a red dashed line.By thorough investigation of this plot, it can be observed that the coherence between waveforms is less guaranteed when specific force is changed quickly.Therefore, the LQR-based MCA can better track and follow the input signals during uniform conditions.The NOC-based technique [17] has a poor COCC of 0.8 and more RMSD of 0.455 than the LQR-based method, which has an improved COCC of 0.88 and minimized RMSD to 0.218.This shows that the LQR method provides much better results in the form of generating realistic motion cues.RMSD should be decreased to zero to achieve realistic motion cues.Fig 8. shows the comparison between normalized lateral specific force inputs and generated motion cues using the LQR-based optimal MCA.Like the previous case of longitudinal specific force, it is observed that the coherence between waveforms is less guaranteed when the specific force is changed quickly.Therefore, LQR-based MCA can better track and follow the input signals during uniform conditions.The NOC-based technique [17] has a lower COCC of 0.83 and RMSD of 0.196 than the LQR-based method, which has an improved COCC of 0.86 and increased RMSD to 0.221.This shows that the LQR method provides much better results in the form of generating realistic motion cues.Greater RMSD shows the doubtful performance of the algorithm during sudden changes of input like braking.To achieve realistic motion, a more precise MCA must be developed that minimizes RMSD to zero.
Fig 9 shows the comparison between yaw velocity and generated velocity cues based on the LQR-based optimal MCA.In contrast to the numerical optimal MCA, this method has an improved COCC of 0.89 and an angle of 3.22.The RMSD was minimized to 0.219.In the case of NOC, these values are 0.81, 2.23, and 0.417, respectively.This shows that LQR-based MCA can decrease sensation error between virtual surroundings and real drivers.The improvement is significant compared with the previous MCA, but it still needs to be improved.
We now compare two popular motion cueing methods based on LQR and GA.The LQR method discussed earlier in this section shows some improvements in RMSD and COCC when compared with the NOC strategy, but it still has undesirable RMSD.All methods are based on the mathematics derived in the motion cueing section, which employs human perception and the human body model to achieve desired results.
These methodologies improve sensation signals by minimizing RMSD and improving the COCC.Using more simulator workspace results in increased sensation error, but an optimal solution may result in the optimized workspace with RMSD as minimum as possible, which needs to have more constraints to be considered with relative simplicity in each constraint.
For GA implementation, some parameters have to be chosen, for example, population size, crossover rate, and mutation rate.By choosing these parameters, the selection is made between the degree of freedom of solution and speed of convergence of the GA.Thus, proper selection of crossover and mutation rate should be ensured to obtain the best results [16].Depending on different applications, these parameters may vary as they are chosen by trial and error.In the case of a large mutation rate, the individual chromosomes may jump to the closer solution, whereas in the case of a small mutation rate, they may stick in local minima.
The recommended settings of these parameters are included in [22].GA consists of the settings of crossover rates and mutation rates, which are satisfactory and meet the demands of the experimental study.Evolutionary operators are assigned to individual chromosomes that constitute the population.It should be noted that the search capacity of the GA is limited because of the small number of chromosomes in the population.If the number of chromosomes increases, the computational time also increases, not suitable for the proposed application.Therefore, an optimal value for population size should be selected to avoid the slow performance of the algorithm.From [7], the population size of fifty chromosomes is appropriate for acceptable results.If the number of chromosomes increases, it will not improve the solution considerably; however, it will significantly decrease computational speed and increase the time of convergence of the solution.
Alternatively, chromosomes are initialized based on the values of the LQR method, as discussed in the previous section.According to the literature [7], the parameter settings for GA parameters are as follows: the crossover rate is 0.8, the mutation rate is 0.2, and half of the size of the chromosomes for individual regeneration.
As the population size tends to increase, the results become more satisfactory in terms of sensation error, but the computational speed tends to decrease accordingly.According to a convergence test on 20, 50, and 100 chromosomes in [16], the best fitness function value, human sensation error, and displacement are obtained using a population size of 50 chromosomes instead of 20.No fruitful results are achieved by increasing the population size from 50 chromosomes.Finally, GA has been successful in decreasing human sensation error and displacement and improving the COCC.
It must be kept in mind that the speed of computation for a population size of 100 chromosomes decreases significantly, resulting in an increased computational burden.When the population size is selected greater than 50, the speed of computation becomes too slow without significant improvement in fitness function value.
Computational time per iteration in the case of 20 chromosomes is 13.90s, whereas in the case of 50 chromosomes is 27.72 s and in the case of 100 chromosomes, it is 52.06 s [16].GAbased optimal MCA stops searching at 455th iteration in case of 50 chromosomes.The density of genetic problems is based on several constraints, generations, and iterations per second, tends to increase with the increase in population size, and is less significantly affected using a population size of 50 compared to a population size of 100.According to Fig 10, the GA-based optimal MCA shows better performance than the LQR method because the GA can track input-specific force with high accuracy.In addition, false cues are removed using the GA.
A comparison of RMSD and COCC between various algorithms is shown in Tables 1 and 2, respectively.The RMSD has been compared with NOC [17], LQR [16], GA [33], OPT [35], and MPC [36].From the Table 1 is can be seen that the NOC behaves imperfectly in case of longitudinal acceleration and yaw velocity, whereas lateral acceleration is inefficiently in case of MPC.This is because the NOC is more related to conventional techniques and MPC cannot better predict the model while in lateral direction.The proposed technique is well suited in all cases as it employs the extended sensation error.Similarly, while we see the coefficient of cross correlation in Table 2 for the sake of comparison between the above techniques, our results are comparable to the GA based approach in all cases of longitudinal & lateral acceleration and also in yaw velocity.This is because our methodology is more closely related to the OPT and GA, benefitting from both.Later we also plot these values of COCC in Fig 10 for the sake of visual comparison with MPC, LQR, GA based techniques.
According to Fig 10, the sensation error between virtual surroundings and the real driver has decreased in the case of the GA method compared to the LQR method.This error can be further decreased and taken below the human threshold using other techniques such as one utilized in the next subsection using H1 control.
A comparison of maximum displacement and maximum angle between various algorithms is shown in Tables 3 and 4, respectively.As shown in Table 3, it is clear that displacement of the simulator in surge direction using GA-based optimal MCA significantly drops to 0.37, which was 0.44 in the case of LQR-based optimal MCA.However, the OPT based approach results least values of both surge and Heave.This can be due to the constraints applied by their methodology.Moreover, the yaw angle from Table 4 in the case of GA-based optimal MCA decreases to 2.04, which was 3.22 in the case of LQR-based optimal MCA.Hence, the GAbased optimal MCA reduces workspace usage by decreasing RMSD and improving COCC.

H1 control-based optimal MCA
The human body model and human perception model have been integrated into the classical washout filter to observe the sensation error, the error between virtual surroundings and real driver.Driver's health comfort and motion sickness are significantly apparent if this error is more than a specified threshold.To minimize this error, hence, to produce realistic motion cues, optimal motion cueing is applied using different optimization methods, like LQR, GA, PSO, and H1.
All previous methods result in some acceptable improvement in error with less workspace usage.H1 control strategies are helpful for their ready applicability to systems involving multiple variables with cross-coupling between channels.However, for an H1 controller, it is difficult to optimize robust performance and stabilization simultaneously.However, many methods, known as extensions of this strategy, exist, such as H1 loop-shaping, which permits the designer to apply classical loop-shaping concepts to the frequency response of the multivariable system and achieve robust performance.
Then favorable robust stabilization is achieved by optimizing the system's response near bandwidth frequency.Here, the Riccati equation is reformulated using optimization and linear matrix inequalities with fewer constraints.show the results of H1 control-based optimal MCA in the case of lateral accelerations, longitudinal acceleration, and yaw velocities, respectively, compared to reference signals along with their cues.
The sensation error between the simulator and driver and the time derivative of the sensational reference signal is input to the H1 control method.Then, the H1 controller generates a supplementary reimbursing signal "Fi," which removes the sensation error between the simulator and driver.The human body has limited ability to enumerate motion applied to it, and the human vestibular system cannot differentiate between small, medium, or large motions in various directions.This motion is divided into logical perceptible ranges.Therefore, in this proposed H1 control method, sensation error values have been divided into five groups: heavily positive, positive, zero, negative, and heavily negative.In this way workspace of the simulator is improved.
In the translational channel of the H1 controller, both inputs are classified as the otolith organ and vestibular organ.The result of the high-pass filter is input to this block and then compared with vehicle input.Then, this error is compared as an input to H1 controller that eliminates error in translational sensational signal.
In the rotational channel of the H1 controller, inputs are provided through the otolith organ and vestibular organ.The result of the low-pass filter is input to this block for tilt coordination and then compared with linear vehicle input to determine the error.Then, this error is then provided as an input to the H1 controller.
As seen in Table 1, in the case of GA-based MCA, the RMSD in the longitudinal direction is 0.072 and the COCC is 0.95.In the case of H1-based optimal MCA, RMSD is significantly reduced to 0.022 and the COCC is significantly improved to 0.98.Moreover, in the case of GA-based optimal MCA, the RMSD in the lateral direction is 0.092 and the COCC in the same direction is 0.96.In the case of the H1-based optimal MCA, RMSD is significantly reduced to 0.036 and the COCC is significantly improved to 0.97.The RMSD in the yaw velocity direction in the case of GA-based optimal MCA is 0. 082 and the COCC in the same direction is 0.98.In the case of the H1-based optimal MCA, RMSD is significantly reduced to 0.013 and the COCC is significantly improved to 0.98.
The RMSD, COCC, maximum displacement, and maximum angle of NOC, LQR, GA, and H1 control-based MCAs are given in Tables 1-4, respectively.
From Table 3, H1-based optimal MCA shows a significant increase in displacement of the simulator to 0.50 in the surge direction.In the case of LQR-based optimal MCA, it is 0.44; in the case of GA-based optimal MCA, it is 0.37.Also, the yaw angle is increased to 4.34 in the case of H-based optimal MCA, 3.22 in the case of LQR-based optimal MCA, and 2.04 in the case of GA-based MCA.Hence, the H1-based optimal MCA improves workspace usage by decreasing RMSD and improving the COCC.This leads to effective implementation of the H1-based optimal routine for motion cueing of driving simulators, considering fewer constraints, and resulting in more workspace usage.
The proposed MCA confirms the production of realistic motion cues by minimizing RMSD.Workspace usage is maximized in two ways: 1. Optimally solving forward kinematics

Using H1 controller
This leads to the effective implementation of optimal routine for motion cueing of driving simulators, considering fewer constraints and resulting in more workspace usage.
In short, all the current optimal MCAs do not utilize a sufficient simulator workspace.Linear data fitting techniques for motion cueing problems are limited by second-order performance indices, such as acceleration values.In contrast, model predictive control involves different systems with different physical constraints.NOC overcomes these limitations by integrating time-varying cost functions adjusted by driver's feedback.Although NOC is an open-loop method, it is still helpful in analyzing performance constraints in loop strategies.Other methods, such as LQR-based optimal control and GA-based optimal control, show better performance in improving the sensation signal.However, they cannot enhance workspace usage due to more constraints involved.The LQR-based method improves the workspace but has more sensation error, and the GA-based method decreases error but with poor workspace usage.

Conclusion
This research proposes a novel optimal MCA for vehicle driving simulators that decreases the sensation error between virtual surroundings and real drivers.The quadratic cost function of the error between vehicle signal and simulator output is minimized using the conventional H1-based optimal control.In addition, usage of the simulator workspace is improved by implementing driver's body model with driver's perception model in optimal MCA.The procedure to design optimal MCA W(s) and the performance of this algorithm through computer simulations have been presented.The simulation results have shown significant improvement both quantitatively and qualitatively.
For future work, other optimization techniques like PSO can be applied for assessment.Moreover, fuzzy logic controller-based optimal MCA can be employed.Like the H1 controller, a fuzzy logic controller can be adopted to enhance motion fidelity.

Fig 3 .
The calculated limb lengths of the Stewart platform are obtained using inverse kinematics, as shown in Fig 4. The method of forward kinematics using the Levenberg-Marquardt method is applied to get values of position and orientation back from values of limb lengths of Fig 4.These values are compared with the previous position and orientation with negligible error.

Fig 13 .
Fig 13.Yaw velocity and resulting cues using the proposed approach.https://doi.org/10.1371/journal.pone.0290705.g013 case of linear motion; A P, Position of moving platform; A a i , Position of A i with respect to O A ; A b i , Position of B i with respect to O A ; B b i , Position of B i with respect to O B ; A i , i th spherical joint at the base platform of Stewart platform; a x , Linear acceleration along the surge; B i , i th spherical joint at moving platform of Stewart platform; d k , Search direction identifier in Gauss-Newton method; dof, Degrees of freedom of robot; e, Error between sensational signal and vehicle signal; f x (s), Force per unit mass along surge; G s , Transfer function for acceleration transition; G T , Transfer function for velocity transition; G Ω , Ratio of vestibular organ to rotational velocity; G ρ , Ratio of otolith organ to the specific force; G BL , Band-limiting filter; G HP , Highpass filter; G LP , Low-pass filter; J, Quadratic cost function; j, No. of joints of a robot; L, Vector of given limb lengths; l i , i th limb length of i th kinematic chain of Stewart platform; n L , No. of links of a robot; O A , Coordinate frame attached to base platform; O B , Coordinate frame attached to moving platform; P i , i th prismatic joint between fixed and moving platform; p x , x-coordinate of position of moving platform; p y , y-coordinate of position of moving platform; p z , z-coordinate of position of moving platform; Q, Weighting matrix for error e; R, Weighting matrix for input signal u s ; t, Vector of position and orientation in screw axis notation; u s (s), Input to driver's body model; u v (s), Input to simulator; v yaw , Yaw velocity; W(s), Transfer function for optimal motion cueing; λ k , Search direction identifier in Levenberg-Marquardt Method; Ω p , Angular velocity or derivative of pitch; Ω ^p(s), Sensation signal in case of angular motion; dof i , Degrees of freedom of i th joint; dof p , Degrees of freedom of passive mechanism; dof s , Degrees of freedom of space; A B R, Orientation of moving platform in rotation matrix notation.